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Abstract 

A Monte Carlo-maximum entropy calculation of the optical conductivity 
of the infinite-dimensional Hubbard model is presented. We show that the 
optical conductivity displays the anomalies found in the cuprate supercon- 
ductors, including a Drude width which grows linearly with temperature, 
a Drude weight which grows linearly with doping, and a temperature and 
doping-dependent mid-IR peak. These anomalies arise as a consequence of the 
dynamical generation of a quasiparticle band at the Fermi energy as T ^ 0, 
and are a generic property of the strongly correlated Hubbard model in all 
dimensions greater than one. 
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I. INTRODUCTION 

The discovery of high-Tc superconductors based on CuO-compound^ has led to a large 
amount of theoretical work about the peculiar properties of these materials. A major ef- 
fort has focussed on the normal-state properties of these compounds. This research was 
largely motivated and substantiated by experiments that revealed striking anomalies in the 
normal-state propertieJi. Most prominent among these are the linear resistivity, a linear 
NMR-relaxation of the Cu-spins, and a Hall angle that behaves ~ over a rather wide 
temperature region. Furthermore, the optical conductivity shows a Drude peak with a 
widthi'S 1/r ~ T, consistent with the linearity of the resistivity, and a Drude weight which 
grows linearly with doping, consistent with the notion of holes acting as the charge carri- 
ers. In addition, there is a pronounced temperature and doping dependent mid-IR peak at 
frequencies above the Drude peak. 

It was argued from the beginning that most of these anomalous properties can be 
explained by two special features appearing simultaneously in these materials: (i) They 
are strongly correlated, i.e. their (effective) local Coulomb interaction is comparable to or 
larger than the characteristic kinetic energy of the relevant carriers; and (ii) they are highly 
anisotropic with the electrons being in principle confined to the CuO-planes characteristic for 
these compounds. Furthermore, it has been argued that the CuO planes can be accurately 
described by a planar single-band Hubbard modeli. The great interest in this class of 
materials has led to a number of new theoretical conjectures, that, although based on the 
assumption of strongly correlated carriers, focussed mainly on the 2D-character of the CuO- 
planesiHli. 

In two earlier publicationsS, we found that several of these anomalies can be understood 
from a Kondo-like effect in the infinite-dimensional Hubbard model. In particular, the 
density of states develops a sharp peak at the Fermi surface as the temperature is lowered. 
The development of this quasiparticle peak coincides with the screening of the effective local 
moments and anomalies in the transport. For example, the resistivity of the model displays a 



2 



distinct linear in T behavior with a slope that increases in inverse proportion to the doping 5, 
consistent with experiment^. The NMR relaxation rate 1/Ti displays a pronounced linear in 
T region, with a slope that increases with doping, also consistent with experimenti. Finally, 
the qualitative features of the Hall resistivity are consistent with experimeniji, including a 
Hall angle that increases quadratically with the temperature. Recently, a sharp Kondo-like 
peak at the Fermi surface has also been seen in the two-dimensional Hubbard mo deli as 
the temperature is lowered. This suggests that these anomalous normal-state properties are 
intrinsic to the Hubbard model, independent of dimensionality! 

In this contribution, we address the anomalous normal-state properties of the Hubbard- 
model optical conductivity. The optical conductivity cr(ct;) is an important probe of the 
excitations of a strongly correlated system. It measures the rate at which electron- hole pairs 
are created by photons of frequency u. In a perfect (translationally invariant) metal, photons 
couple only to electron-hole pairs with vanishing momentum and energy; o'{uj) is proportional 
to a Dirac delta function [cr(u;) = Dperfect^i^^)] with Drude weight Dperfect = 'ne^n/m (we 
set h = 1). Electron-electron correlations modify this picture at zero temperature: The 
charge and spin fluctuations induce a dynamic disorder to the lattice potential which reduces 
the free-carrier Drude weight by the inverse of the quasiparticle renormalization factor Z 
{D = Z^^Dperfect) and transfers the remaining spectral weight to a frequency-dependent 
component of o"(a;) reflecting the incoherent charge and spin fluctuations; the total spectral 
weight is, however, conserved. Finite-temperature effects will broaden the zero-frequency 
delta function into a Lorentzian and can modify both the quasiparticle renormalization and 
the higher-frequency excitations. This simple picture is further modified when restriction 
is made to a single (or finite number) of electronic bands. In this case, the total spectral 
weight can vary as a function of temperature or interaction strength because the projection 
onto a restricted basis set disregards all excitations to electronic bands that are higher in 
energy than those that have been kept in the model. 

If the on-site Coulomb repulsion U is large enough, then the system will be a Mott 
insulator at half-filling. The single-particle density of states N{uj) consists of two symmetric 
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bands, separated by U (called the lower and upper Hubbard bands) with the Fermi level 
lying in the middle of the gap. The renormalization factor Z diverges since there are no 
quasiparticles at the Fermi level. The optical conductivity will consist of a charge-excitation 
peak centered at u; ~ f/, whose width is the order of the bandwidth. As the system is doped 
away from half filling, a quasiparticle resonance appears at the Fermi energy as T ^ 0. The 
weight of the quasiparticle peak increases with doping. In this case, the optical conductivity 
will have a Drude peak (from the "free" quasiparticles at the Fermi energy) and a mid-IR 
band (because of excitations between the lower Hubbard band and the quasiparticle peak) 
in addition to the charge-excitation peak. The Drude weight can be estimated as follows: 
when S holes are doped into the half- filled band, the total electron concentration is {1 — 6), 
and the quasiparticle renormalization factor behaves roughly as Z oc 1/5, so the Drude 
weight is expected to be 

Dorude OC Z^'^n (X 6{1 - 6) , (1) 

which increases linearly for small and large doping, and is peaked near quarter filling {n = 
0.5). In addition, the weight of the mid-IR peak should also behave qualitatively like Eq. (|Tp 
because it involves excitations between the quasiparticle peak and the lower Hubbard band, 
which also should grow as the density of holes in the half-filled band. The "hole-like" nature 
of the charge excitations near half filling results from the strong Coulomb renormalizations 
that create the Mott insulator at exactly half filling. 

In this contribution, we investigate the optical conductivity of the Hubbard model in in- 
finite dimensions. Our methodology is detailed in Section II, where we discuss the formalism 
and calculational techniques. Section HI contains our results, which include a Drude peak 
whose width 1/r grows linearly with temperature, a mid-IR peak that becomes more visible 
at low temperatures and dopings, and a charge-excitation peak that decreases with doping 
and is weakly dependent upon T. Comparison of our results with experiment is presented 
in Section IV, and conclusions are given in Section V. 
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II. METHOD 

Motivated by the observations of Anderson! and Zhang and Rice! we will study the 
single-band Hubbard modefl in d-dimensions, with 

t* 

H = -— ^ E 4^- + E ^i^^ii ' (2) 

as a model for the CuO system. Our notation is the following: is an electron creation 
operator for an electron in a localized state at lattice site i with spin cr; Uia- = cl^Qa is the 
corresponding electron number operator; t* is the rescaled hopping matrix element; and U 
is the Coulomb interaction strength. We choose t* = 1 as a convenient unit of energy. This 
model may be solved in the limit of high dimensions, using the observations of Metzner and 
Vollhardtlll that with increasing coordination number of the underlying lattice the many- 
body renormalizations due to a two-particle interaction like the Hubbard-[/ in (0) become 
purely local. 

The solution of the model may be mapped onto the solution of purely local correlated 
system coupled to an effective (self-consistently determined) bathlil^. The quantum Monte 
Carlo (QMC) scheme based on the work of Hirsch and Fye^Ilil has proven to be the most 
effective method for solving this strongly correlated local problem. Dynamical properties of 
the model are then obtained using numerical analytic continuation (employing maximum- 
entropy techniques)^. This method requires a default model. To obtain the single-particle 
density of states we use the finite-U non-crossing approximation (NCA)@ result for the 
infinite-dimensional Hubbard model as a default model at high temperatures, where it is 
essentially exact0. At low temperatures, where the NCA fails, we use the numerically 
continued result of the next higher temperature as a default model. The posterior probability 
of the final result is employed to determine which default model should be used. Generally, 
we find that the crossover temperature between using the NCA default model and the 
higher-temperature continuation lies at T ^ 2Tq, where Tq is the Kondo-like scale for this 
mo del§. Once the density of states is determined, the self energy S(u;) may then be found 
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by (numerically) inverting the Fadeev function w{z) in the relation 

N{uj) = Re{w[uj + fx-J:{w)]} /^/^. (3) 



With the knowledge of the one-particle self energy, one can calculate transport quantities. 
For example, the optical conductivity can be calculated exactly in the local approximation. 
It is given by the simple bubble onlyil'il, whose evaluation leads to 



— oo k(7 ^ ^ 



oo 



= arc J de^-^^^-^^^±^ j dyp{y)A{y,e)A{y,e + u) , (4) 

— oo — oo 

where V is the lattice volume, a = e^/ (Vd) and defines the unit of the conductivity, the 
spectral weight satisfies A(ek,io) = —-^^^ [G{k,u)], and the noninteracting density of states 
is p{y) = exp(-y2)/^. 

As the temperature is lowered, the Hubbard model in infinite-d is found to always be a 
Fermi liquid^S, except for the region of phase space where it is magneticS. A Fermi liquid 
is defined by a self-energy that has the following structure: 

ReE(w + iO+) = ReE(O) - Z) + 0(cj2) , 

ImS(a; + iO+) = -r + 0(cu2) , (5) 

with r oc for temperatures T <^ Tq the characteristic Fermi temperature. [The Fermi 
temperature Tq decreases to zero as half filling is approached^, and the Fermi-liquid-theory 
form of Eq. still holds for moderate temperatures, with the only change being F oc T for 
T > Tq.] The spectral weight, then assumes the form 

1 r 

A{y,uj) = - , . y , - + Ajr,c{y,^) , (6) 

TT 1 ^ + [uZ + ep - yy 

with the Fermi level defined hj tp = p — ReS(O) and Ainciy^^) denoting the (rather struc- 
tureless) incoherent contributions to the spectral function. The spectral function includes 
a delta function at zero temperature [Aixj^oj) 5{ijjZ + ep — y) + A/„c(l/,^^)] because the 
broadening F vanishes in that limit. 
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If the Lorentzian form for the Fermi-hquid-theory spectral function [Eq. (^] is substi- 
tuted into the expression for the optical conductivity found in Eq. and the temperature 
satisfies T ^ Tq ^ t*, then the optical conductivity becomes 

a vr 1 + cu^r^ a 
with Drude weight D = Z^^np^ep), relaxation time r = 1/{2T) oc and <7inc{^^) 

containing the contributions from the incoherent pieces of the spectrum. 

The noninteracting Drude weight for the single-band model {Z = 1) satisfies -D„on = 
np^ep)- As the interaction U is turned on in a single-band model, the total integrated 
spectral weight is not conserved but becomes [/-dependent^ because spectral weight that 
would be shifted to higher bands is "lost" in any single-band model. The total spectral 
weight satisfies 

roo 

2 / a{uj)duj = -vre^ (T,) . (8) 
Jo 

where (— T^) is the kinetic energy per site, divided by the number of lattice dimensions. It 
should be stressed that the modification of the sum rule in Eq. (^) will produce some system- 
atic modifications to the behavior of the optical conductivity as a function of temperature, 
interaction strength, and doping. These systematic effects must be kept in mind when one 
is comparing the results of a single-band calculation to experiment. 

The Drude weight D for the interacting system, may also be determined (at T = 0) 
by extrapolation of the Matsubara-frequency current- current correlation function using the 
method proposed by Scalapino et al.0 {D measures the "free" quasiparticles in the system). 
This method determines the Drude weight of a metal by examining the asymptotic form 
of the current-current susceptibility in the x-direction, A^xi^, iuJn)- More specifically, D is 
given by 



D = lim TT 

T->0 



e2(-T,)-A,.,.(q = 0,2mT) , (9) 



where the limit T ^ is taken after first setting the momentum transfer to zero (q = 0). 
Note that this latter method of determining the Drude weight is a much better defined 



procedure than trying to fit the optical conductivity to the generic form of Eq. (0) because 
of the uncertainty left in trying to fit ajnc{uj). 



III. RESULTS 

We present here our results for the optical conductivity of the single-band Hubbard model 

in infinite dimensions with U = At*. Figure |l](a) shows the optical conductivity obtained 

from Eq. when 6 = 0.068 for a variety of temperatures. One finds the Drude peak at 

uj = developing with decreasing temperature. In addition, a shoulder develops adjacent to 

the Drude peak at ~ 1 which is strongly temperature dependent and clearly visible only 

for the lowest temperatures. The last feature in (t{uj) is a roughly temperature-independent 

peak a.t UJ ^ U. In order to compare our results to experiment, the three features in (t{uj) are 

fit to a Lorentzian plus (asymmetric) harmonic-oscillator forms for the higher-energy peaks 

crjuj) D T Cmir (^'^'^mir Cg uj'^Tq , . 

a ^ TT 1 + co'r^ + n uj^Tl,,^, + (uj^ - + tt u^Vl + (o.^ - u^lY ' ^ ^ 

with r the relaxation time of the quasiparticles, and the constants Cmir, ^mir, and Tmir 

{Cc, ujc, and r^) the weight, center, and width, respectively of the mid-IR (charge-transfer) 

peak. The Drude width 1/r obtained from this fitting procedure is plotted in the inset to 

Fig. [5(b). Note that for temperatures on the order of To, it is well approximated by a 

straight line (the line is a guide to the eye), while for T <^Tq the Drude width must change 

its behavior to r oc according to the general properties of a Fermi liquid. By subtracting 

off the fit Drude portion from the optical conductivity, we were able to isolate the mid-IR 

portion, as shown in Fig. |l](b). Note that the "double-peak" structure emerging in the 

mid-IR peak at the lowest temperature is an artifact of the fitting procedure which is not 

perfect in extracting the parameters for the Drude peak. The mid-IR peak is temperature 

dependent, growing in size and moving to slightly lower frequencies as the temperature is 

lowered. 

Comparison of the Drude weight D determined by the fitting procedure in Eq. (|TUp 
and the independent method of calculation in Eq. (^) produces only qualitative agreement. 
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We have tried more sophisticated fitting routines in which we fit the Fermi-hquid-theory 
parameters {Z, T, and e^) in the spectral function of Eq. and then determine the Drude 
contribution to the optical conductivity by employing the full expression in Eq. (|^), and the 
fit improves, as does the comparison of the Drude weight, but the quantitative agreement 
still has a systematic error on the order of 10 — 20% (where the fitting procedure always 
overestimates the Drude weight), because of the simple form chosen for the mid-IR peak. 
Our conclusion is that the method of Ref. ^ is superior to any ad hoc fitting procedure 
in determining the low-temperature Drude weight. Note, however, that the charge-transfer 
peak is uniquely determined by the above fitting procedure (in the sense that the results are 
independent of what fitting procedure is employed), indicating that the harmonic-oscillator 
form is a reasonable approximation to that peak. 

The existence of the Drude and charge-excitation peaks have been reported previouslyS. 
These results were obtained with the NCA and were thus restricted to temperatures T ^ 2To. 
The mid-IR bump, however, becomes clearly visible only for temperatures below this scale. 
We can presently resolve this feature because of the refined numerics employed in the QMC 
and the numerical analytic continuation at lower temperatures. 

In Fig. the optical conductivity is plotted as a function of doping for fixed temperature 
/3 = 43.2. The different dopings are best identified by their decreasing charge-transfer peaks 
(S = 0.068, 0.0928, 0.1358, 0.1878, 0.2455, 0.3, 0.35, 0.4, and 0.45). The solid lines indicate 
the electron concentrations where the optical conductivity increases with doping in the lower- 
Hubbard-band region of u < 2 {6 < 0.25) and the dotted lines are where (t{uj) decreases with 
doping in the same region {6 > 0.25). One can see that the mid-IR peak is strongly doping 
dependent, being most distinct from the Drude peak at low dopings, and merging with it 
as the doping and the width of the Drude peak increases. In the inset, we show that the 
Drude weight D increases linearly with 6. This latter result shows that the "free" carriers in 
a doped Mott insulator are holes in the half-filled band. As the doping is increased further, 
the Drude weight eventually saturates, and then decreases with doping, as the character of 
the "free" carriers changes from being hole-like to being electron-like. 
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The optical conductivity clearly displays an isobestic point, where cr(a;) is independent 
of doping {ujjB ~ 2). The isobestic point marks the boundary between the regions where 
the weight of the optical conductivity increases as a function of doping < ujib), and 
the regions where the weight decreases with doping {u > ujjb)- However, it is incorrect to 
assume that all of the spectral weight that lies below the isobestic point was transferred 
from above. An effective carrier number Neff{uj), defined by the integral of the optical 
conductivity 



is plotted in Figure ^ (the normalization is chosen here to give the number of carriers for a 
model with a complete set of basis states in order to make contact with experiment). The 
integrated optical conductivity increases rapidly with doping for small values of 6, but then 
increases more slowly, eventually saturating, and decreasing with doping (at 6 ~ 0.4). Once 
again, solid lines denote the electron concentrations where Neff{uj) increases with doping 
{6 < 0.4) and the dotted line is where N^ffioj) decreases with doping (5 > 0.4). Note that 
there is no frequency where the integrated spectral weight is independent of doping, which 
would be required for the above scenario (of the Drude and mid-IR weight coming entirely 
from the charge-transfer peak) to occur. 

In addition, we can examine how the sum rule for the total spectral weight, and for 
it's component pieces (Drude, mid-IR, and charge-transfer) evolve as a function of doping 
at fixed low temperature (in fact, the following analysis is an approximation to the zero- 
temperature behavior of the different components of the spectral weight). This is shown in 
Figure ^ The total spectral weight is found from Eq. (^ (which agrees with the integral of 
the optical conductivity to within 3%), the charge-transfer weight from the fitting procedure 
in Eq. (|10|), the Drude weight from Eq. (^), and the weight of the mid-IR peak is then 
determined by subtraction (the charge-transfer weight is neglected for n < 0.5). Note how 
the mid-IR peak initially increases linearly with doping until n ~ 0.75 where it saturates and 
then decreases as the system becomes less correlated, and how the system changes from a 




(11) 
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hole-like metal to an electron-like metal at about quarter filling. In addition, one can easily 
see that the Drude and mid-IR weights both grow quite rapidly as the system is doped away 
from half filling (the initial growth rate of the mid-IR weight is about a factor of 2 faster than 
that of the Drude weight). This fast growth of the Drude and mid-IR weights explains why 
the optical conductivity for u < ujb increases so rapidly with doping as shown in Fig. |^. 
These results are qualitatively similar to those found in the strong-coupling limit of the 
Hubbard model in one dimension^, indicating that the optical properties of the Hubbard 
model are not strongly dependent upon the dimensionality (with the exception that there 
is no mid-IR peak in one dimension). 

One can also inquire into the temperature dependence of the Drude, mid-IR, and charge- 
transfer peaks at fixed doping. The analysis, in this case, relies on trying to fit the optical 
conductivity to the form in Eq. |T0|, as a function of temperature. Such a procedure is 
ill-defined, because the Drude and mid-IR peaks merge at high temperatures, and it is 
difficult to separate them into their component pieces, without explicit knowledge of the 
proper fitting forms. However, the qualitative physics can be discussed: We find that as 
the temperature is raised, that spectral weight is lost from the system (i. e., the magnitude 
of the expectation value of the kinetic-energy operator decreases as T increases), that the 
Drude and charge-transfer weights slightly increase, and that the mid-IR weight is sharply 
reduced. In addition, spectral weight rapidly shifts to higher frequency, since the Drude 
width (which is easy to determine by the fitting procedure) increases linearly with T for 
temperatures larger than Tq. 

The origin of all of these features in the optical conductivity can be understood by 
studying the single-particle density of states, plotted as a function of both temperature (a) 
and doping (b) in Fig. ^. Clearly the Drude peak results from the quasiparticle band which 
develops at low temperatures. We have shown previously that while this peak is developing 
[i. e. in the region T > Tq), the scattering rate (as measured by the resistivity) increases in 
proportion to the temperature^. Thus, the Drude width also increases in proportion to the 
temperature. We attribute the mid-IR peak to excitations from the lower Hubbard band to 
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the quasiparticle band at the Fermi energy. Note that the quasiparticle resonance is only 
sharp and distinct from the lower Hubbard band at low T and small 6. Thus, the mid-IR 
peak will only be visible in the optical conductivity at low T and small 6 as well. On the 
other hand, the upper peak at ^ 4 is due to the charge excitations from both the lower 
Hubbard band and the quasiparticle band to the upper Hubbard band. Since the upper 
Hubbard band is distinct from the quasiparticle weight at the Fermi surface, and since the 
weight in the upper Hubbard band does not change as significantly with doping, the upper 
peak in cr(to') has a weaker temperature and doping dependence, but is expected to disappear 
as the system becomes uncorrelated in the low-density limit. 



Many experimental measurements have been made of the optical conductivity in doped 
cuprates. Kramers-Kronig analysis is used to determine cr(ci;) from reflectivity measurements 



ini La2-xSr^Cu04 (LSCO), iniH YBaaCuaOe+x (YBCO), inill BiaSrsCaCuzO^ (BSCCO), 



is also employed to measure the optical conductivity inEl YBCO and in hotim LSCO and 
NCO. It is interesting to compare the theoretical calculation of the optical conductivity with 
these experimental results. 

The experiments yield six trends for the cuprate superconductors: 1) The mid-IR peak 
maximum moves to lower frequency, and merges with the Drude peak as the doping increases; 
its spectral weight grows very rapidly with doping near half filling; 2) At a fixed value of 
the doping, spectral weight rapidly moves to lower frequency as T — 0, but the total weight 
in the Drude and mid-IR peaks remains approximately constant; the width of the Drude 
peak decreases linearly with T; 3) The insulating (or undoped) phase has a charge-transfer 
gap; when doped, the optical conductivity initially increases within the gap region, but 
eventually saturates and then decreases with doping; 4) There is an isobestic point or nearly 
isobestic behavior (in that the optical conductivity is nearly independent of doping) at a 



IV. COMPARISON TO EXPERIMENT 



and inSi Nd2Cu04-j, (NCO) 
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frequency that is approximately one half of the charge-transfer gap; 5) The effective charge 
has a constant value with respect to doping near the charge-transfer edge; and 6) more than 
one peak is observed within the mid-IR region. 

Most of these trends are reproduced by our theoretical model. In particular, Fig. |^ 
illustrates how the mid-IR peak moves to lower frequency and joins the Drude peak (1), how 
the optical conductivity initially increases with doping at low frequency, but then saturates 
and decreases (3), and how there is an isobestic point (4). Fig. || shows how the spectral 
weight is transferred to lower frequencies as the temperature is lowered, but we find that 
the total (Drude plus mid-IR) spectral weight does not remain constant as in (2) because 
the temperature dependence of the expectation value < > produces some temperature 
dependence to the total Drude plus mid-IR spectral weight (it actually decreases as T 
increases). Fig. |^ does not display the trend of point (5), possibly because the restriction to 
a single-band model reduces the optical conductivity at higher frequencies to such an extent 
that the effective charge must depend on the doping level for the theoretical model. We do 
not see the multiple peaks in the mid-IR region of point (6). 

These experimental features are usually attributed^ to phonons or impurities present 
in the system, but judging from our results, the low-energy feature may also be due to 
excitations from the lower Hubbard band to a dynamically generated quasi-particle band at 
the Fermi energy (which is generated by Kondo-like screening of the moments). Naturally, 
the single-band model fits the experminental data much better at lower frequencies, where 
the single-band approximation is relevant, but fails in reproducing some of the higher-energy 
trends found in the cuprates, because of the neglected bands. A comprehensive theory should 
include both the effects of strong electron correlation (which reproduce the insulator at half 
filling, and give a hole-like character to the charge excitations near half filling) with the 
effect of electron-phonon coupling (which are necessary to explain similar mid-IR features 
in nonsuperconducting perovskites) . 
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V. CONCLUSION 



A quantum Monte Carlo and maximum entropy calculation of the optical conductivity of 
the infinite-dimensional Hubbard model has been presented. The Mott-insulating character 
of the ground state at half filling drives many anomalous behaviors in the normal state near 
half filling that are similar to those observed in the cuprate superconductors. In particu- 
lar the system is always a Fermi liquid away from half filling, but the Fermi temperature 
vanishes, and the quasiparticle renormalization factor diverges as half filling is approached. 
As a result, the free carriers in the system initially have a hole-like character (that changes 
to an electron- like character at approximately quarter filling). The Drude width for these 
carriers grows linearly with temperature for temperatures above To, the Drude weight grows 
linearly with doping, and there is a doping and temperature dependent mid-IR peak. These 
anomalies arise naturally from the presence of a strongly temperature-dependent quasiparti- 
cle peak, whose origin is a Kondo-like screening of the magnetic moments, and which appears 
to occur in the Hubbard model for all dimensions greater than 1. 

The anomalous features in the experimentally measured optical conductivity for the 
cuprates are usually attributed to either polarons or impurities. However, any purely po- 
laronic theory has difficulty in explaining the magnetic insulating character of the ground 
state at half filling. The Hubbard model naturally describes such an insulating state, and 
appears to also describe many of the anomalous features present in the experimental data. 
In light of this fact, it is worthwhile to try to incorporate both the effects of strong electron 
correlation, and the electron-phonon interaction into a comprehensive theory for the normal 
state of the cuprate materials. Work along these lines is in progress. 
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FIG. 1. (a) Optical conductivity vs. u for various temperatures when U = 4 and S = 0.068 in 
units of a = e^ /iVd). Note that at low temperatures, when the Kondo peak becomes pronounced in 
the DOS (see Fig. a mid-IR feature begins to appear in a{uj). The mid-IR feature is isolated 
in (b) by fitting the low-frequency data with Eq. and subtracting off the Drude part. As the 
temperature is lowered, the mid-IR peak becomes more pronounced and shifts to lower temperatures. 
Note that the double-peak structure in the mid-IR peak is most likely due to the inaccurate fitting 
form for aMini^)- shown in the inset, the width (1/t) of the Drude peak is found to increase 
roughly linearly with T. 
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FIG. 2. The doping dependence of the optical conductivity when {7 = 4 and (3 = 43.2. Note 

that for larger 5 the mid-IR and Drude peak begin to merge, so that the latter is less distinct. The 
different dopings are identified by their decreasing charge-transfer peaks, respectively (5 = 0.068, 
0.0928, 0.1358, 0.1878, 0.2455, 0.3, 0.35, O.4, and 0.^5). The solid lines correspond to the case 
where the optical conductivity increases with doping in the lower-Hubbard-band region of to < 2 
{5 < 0.25) and the dotted lines correspond to the case where (t{uj) decreases with doping (5 > 0.25). 
The inset shows the evolution of the Drude weight D as a function of doping, which is computed 
using Eq. ^ at (3 = 20. For 6 ^ 0.4, the Drude weight increases linearly with doping. 
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FIG. 3. Effective charge N(,ff{uj) of the Hubbard model when U = 4 and (3 = 43.2 as a function 

of doping. The same electron concentrations are plotted as in Figure [J. The solid lines denote 
doping levels where the effective charge increases with doping, and the dotted line corresponds to 
the case where the effective charge decreases with doping (5 > OA). 
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FIG. 4. Total spectral weight and its individual components for the optical conductivity of the 

Hubbard model as a function of doping at U = A and (3 = 43.2. The total spectral weight (solid 
dots) is broken up into its component pieces (Drude [open triangles], mid-IR [open squares], and 
charge-transfer [solid triangles]). Note how the charge-transfer peak is important only as the elec- 
tron concentration nears half filling, and how the mid-IR peak disappears as the electron correlations 
become smaller. 
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FIG. 5. One-particle DOS of the Hubbard model. In ( a) the DOS is plotted for different dopings 
6 = 1 — n at 13 = 43.2. In addition to the lower and upper Hubbard bands, which are fairly doping 
independent, a resonance occurs at the chemical potential that becomes broader with increasing S 
and finally merges with the lower Hubbard band. In (b ), the temperature evolution of the density of 
states is plotted for U = A and 6 = 0.188. As the temperature is lowered, a sharp peak, distinct from 
the lower Hubbard band, develops at the Fermi surface. As shown in the inset, the development of 
a sharp peak at the Fermi surface is correlated with the reduction of the screened local moment and 
hence is associated with resonant Kondo screening of the spins. 



24 



